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Cs) , We review lattice QCD results for glueballs (including a discussion of 

£3 \ mixing with scalar mesons), hybrid mesons and exotic mesons (such as 

B S B S molecules). 
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1 Introduction 
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The most systematic approach to non-perturbative QCD is via lattice tech- 
niques. Lattice QCD needs as input the quark masses and an overall scale 
(conventionally given by Aqcd)- Then any Green function can be evaluated 
by taking an average of suitable combinations of the lattice fields in the vacuum 
samples. This allows masses to be studied easily and matrix elements (partic- 
ularly those of weak or electromagnetic currents) can be extracted straightfor- 
wardly. Unlike experiment, lattice QCD can vary the quark masses and can 
also explore different boundary conditions and sources. This allows a wide 
Q-i! range of studies which can be used to diagnose the health of phenomenological 

J|h \ models as well as casting light on experimental data. 

One limitation of the lattice approach to QCD is in exploring hadronic 
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decays because the lattice, using Euclidean time, has no concept of asymp- 
totic states. One feasible strategy is to evaluate the mixing between states 
of the same energy - so giving some information on on-shell hadronic decay 
amplitudes. 

For comparison with models and for ease of computation, the special case 
of infinitely heavy sea quarks (namely neglect of quark effects in the vacuum: 
the quenched approximation) is often used. We shall also present results from 
including sea quark effects - usually two flavours of degenerate sea quark of 
mass equivalent to strange quarks or heavier. 

The quark model gives a good overall description of hadronic spectra. Here 
I will discuss lattice results for states which go beyond the quark model: glue- 
balls, exotic mesons and hybrid mesons. 



2 Glueballs 



Glueballs are defined to be hadronic states made primarily from gluons. The 
full non-perturbative gluonic interaction is included in quenched QCD. In the 
quenched approximation, there is no mixing between such glueballs and quark 
- antiquark mesons. A study of the glueball spectrum in quenched QCD is 
thus of great value. This will allow experimental searches to be guided as well 
as providing calibration for models of glueballs. A non-zero glueball mass in 
quenched QCD is the "mass-gap" of QCD. To prove this rigourously is one 
of the major challenges of our times. Here we will explore the situation using 
computational techniques. 

In principle, lattice QCD can study the meson spectrum as the sea quark 
mass is decreased towards experimental values. This will allow the unambigu- 
ous glueball states in the quenched approximation to be tracked as the sea 
quark effects are increased. It may indeed turn out that no meson in the phys- 
ical spectrum is primarily a glueball - all states are mixtures of glue, qq, qqqq, 
etc. We shall later discuss lattice results on the mixing of glueballs and scalar 
mesons (ie qq states). 

In lattice studies, dimensionless ratios of quantities are obtained. To ex- 
plore the glueball masses, it is appropriate to combine them with another very 
accurately measured quantity to have a dimensionless observable. Since the 
potential between static quarks is very accurately measured from the lattice, 
it is now conventional to use r for this comparison. Here r is implicitly 
defined by r 2 dV(r)/dr = 1.65 at r = r$ where V(r) is the potential energy 
between static quarks which is easy to determine accurately on the lattice. 
Conventionally tq w 0.5 fm. 

Theoretical analysis indicates that for Wilson's discretisation of the gauge 
fields in the quenched approximation, the dimensionless ratio mro will differ 
from the continuum limit value by corrections of order a 2 . Thus in fig. 1 the 
mass of the J PC =0 ++ glueball is plotted versus the lattice spacing a 2 . The 
straight line then shows the continuum limit obtained by extrapolating to 
a = 0. As can be seen, there is essentially no need for data at even smaller 
a-values to further fix the continuum value. The value shown corresponds 
to m(0 ++ )r = 4.33(5). Since several lattice groups 0, |3], f| have mea- 
sured these quantities, it is reassuring to see that the purely lattice observ- 
ables are in excellent agreement. The publicised difference of quoted m(0 ++ ) 
from UKQCD [|3| and GF11 Q comes entirely from relating quenched lattice 
measurements to values in GeV. 

In the quenched approximation, different hadronic observables differ from 
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Figure 1: The value of mass of the J PC = ++ glueball state from quenched 
data (N F = 0)[TJ, ^ |3|, f| in units of r where r « 0.5 fm. The straight 
line shows a fit describing the approach to the continuum limit as a — > 0. 
Results P, |20L El with N F = 2 flavours of sea quarks are also shown. 



experiment by factors of up to 10%. Thus using one quantity or another to 
set the scale, gives an overall systematic error. Here I choose to set the scale 
by taking the conventional value of the string tension, y/a = 0.44 GeV, which 
then corresponds to r^ 1 = 373 MeV. An overall systematic error of 10% is then 
to be included to any extracted mass. This yields m(0 ++ ) = 1611(30) (160) 
MeV where the second error is the systematic scale error. Note that this is 
the glueball mass in the quenched approximation - in the real world significant 
mixing with qq states could modify this value substantially. 

In the Wilson approach, the next lightest glueballs are |2], |3| the tensor 
m (2++)r = 6.0(6) (resulting in m(2 ++ ) = 2232(220) (220) MeV) and the pseu- 
doscalar m(2 ++ )r = 6.0(1.0). Although the Wilson discretisation provides a 
definitive study of the lightest (0 ++ ) glueball in the continuum limit, other 
methods are competitive for the determination of the mass of heavier glue- 
balls. Namely, using an improved gauge discretisation which has even smaller 
discretisation errors than the a 2 dependence of the Wilson discretisation, so al- 
lowing a relatively coarse lattice spacing a to be used. To extract mass values, 
one has to explore the time dependence of correlators and for this reason, it is 
optimum to use a relatively small time lattice spacing. Thus an asymmetric 
lattice spacing is most appropriate. The results [0] are shown in fig. 2 and 



for low lying states are that m(0 ++ )r = 4.21(H)(4), m(2 ++ )r = 5.85(2)(6), 
m(0 _+ ) = 6.33(7)(6) and m(l + -)r = 7.18(4)(7). Another recent study § has 
used an improved discretisation based on the perfect action approach (without 
a space-time asymmetry) and obtains results consistent with earlier work. 

One signal of great interest would be a glueball with J PC not allowed for 
qq - a spin-exotic glueball or oddball - since it would not mix with qq states. 
These states are found |2], 0] to be high lying: considerably above 2m(0 ++ ). 
Thus they are likely to be in a region very difficult to access unambiguously 
by experiment. 

Within the quenched approximation, the glueball states are unmixed with 
qq, qqqq, etc. Furthermore, the qq states have degenerate flavour singlet and 
non-singlet states in the quenched approximation. Once quark loops are al- 
lowed in the vacuum, for the favour-singlet states of any given J PC , there will 
be mixing between the ss state, the uu + dd state and the glueball. One way to 
explore this is to measure directly the scalar mass eigenstates in a study with 
Nf = 2 flavours of sea-quark. Most studies show no significant change of the 
glueball spectrum as dynamical quark effects are added - but the sea quark 
masses used are still rather large || |20| . A recent study |§, however, does 
find evidence for a reduced mass, albeit with a rather large lattice spacing, see 
fig. 1. This effect could be due to mixing of scalar mesons and glueballs, as we 
discuss below, or might just be a sign of an enhanced order a 2 correction at 
the relatively large lattice spacing used. 

Let us now discus the mixing of the scalar glueball and scalar mesons. The 
mass spectrum of qq states has been determined on a quenched lattice and 
the scalar mesons are found to lie somewhat lighter than the tensor states II . 
These 2 ++ mesons are experimentally almost unmixed and so will be quite 
close to the quenched mass determination. This suggests that the quenched 
scalar masses from the lattice are at around 1.2 GeV and 1.5 GeV (for nn and 
ss respectively). An independent study [JlO|, |ll] suggests that the scalar ss 
state is about 200 MeV lighter than the glueball which is a broadly compatible 
conclusion. Thus the glueball, at around 1.6 GeV, lies heavier than the lightest 
qq scalar states. This information can then be combined with mixing strengths 
to give the resulting scalar spectrum. 

It is possible to measure the mixing strength on a quenched lattice even 
though no mixing actually occurs. On a rather coarse lattice (a" 1 ~ 1.2 GeV), 



two groups have attempted this |TT], g. Their results expressed as the mixing 
for two degenerate quarks of mass around the strange quark mass are similar, 
namely E « 0.36 GeV [M and 0.44 GeV @. From this evaluation of the 



mixing strength, one can use a mass matrix to estimate the mass shift induced 
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Figure 2: The continuum glueball spectrum | 



in the glueball and scalar meson. The relevant mass matrix is (in a glueball, 
qq basis in GeV units): 

/ 1.1 0.4 \ 

^ 0.4 1.6 ) 

which would give a downward shift of the glueball mass by 20%. This is in 
qualitative agreement with our direct determination with Nf = 2 flavours of 
sea-quark that the lightest scalar mass is reduced significantly at this lattice 
spacing as shown in fig. 1. 

Note that at this coarse lattice spacing the quenched glueball mass is re- 
duced (see fig. 1) below the canonical value of 1.6 GeV. Thus a study at 
smaller lattice spacing is needed. An exploratory attempt to extrapolate to 
the continuum (jTIJ] gave a very small mixing of 86(64) MeV, while the other 
determination || uses clover improvement so order a effects in the extrapola- 
tion to the continuum are suppressed and one would not expect a significant 
decrease in going to the continuum limit. What this discussion shows is that 
precision studies of the mixing on a lattice have not yet been achieved. 

As well as this mixing of the glueball with qq states, there will be mixing 
with qqqq states which will be responsible for the hadronic decays. A first 
attempt to study this |12| yields an estimated width for decay to two pseu- 
doscalar mesons from the scalar glueball of order 100 MeV. A more realistic 
study would involve taking account of mixing with the nn and ss scalar mesons 
as well. 



3 Exotic states 



By exotic state we mean any state which is not dominantly a qq or qqq state. 
Such examples have been known for a long time: the deuteron is a proton- 
neutron molecule for example. It is very weakly bound (2 MeV) and is quite 
extended. Similar molecular states involving two mesons have been conjec- 
tured. 

One case which is relatively easy to study is the BB system, idealised 
as two static quarks and two light quarks. Then a potential as a function 
of the separation R between the static quarks can be determined. Because 
the static quark spin is irrelevant, the states can be classified by the light 
quark spin and isospin. Lattice results [|13|] (using a light quark mass close 



to strange) have been obtained for the potential energy for I q = 0, 1 and 
S q = 0,1. For very heavy quarks, a potential below 2Mb will imply binding of 
the BB molecules with these quantum numbers and L = 0. For the physically 
relevant case of b quarks of around 5 GeV, the kinetic energy will not be 
negligible and the binding energy of the BB molecular states is less clear cut. 
One way to estimate the kinetic energy for the BB case with reduced mass 
circa 2.5 GeV is to use analytic approximations to the potentials found. For 
example the I q , S q =(0,0) case shows a deep binding at R = which can be 
approximated as a Coulomb potential of — 0.1/ R in GeV units. This will give 
a di-meson binding energy of only 10 MeV. For the other interesting case, 
(I q , S q )=(0,l), a harmonic oscillator potential in the radial coordinate of form 
— 0.04[1 — (r — 3) 2 /4] in GeV units leads to a kinetic energy which completely 
cancels the potential energy minimum, leaving zero binding. This harmonic 
oscillator approximation lies above the estimate of the potential, so again we 
expect weak binding of the di-meson system. 

Because of these very small values for the di-meson binding energies, one 
needs to retain corrections to the heavy quark approximation to make more 
definite predictions, since these corrections are known to be of magnitude 46 
MeV from the B, B* splitting. It will also be necessary to extrapolate the 
light quark mass from strange to the lighter u, d values to make more definite 
predictions about the binding of BB molecules. 

Models for the binding of two B mesons involve, as in the case of the 
deuteron, pion exchange. The lattice study |l3j is able to make a quantita- 
tive comparison of lattice pion exchange with the data described above and 
excellent agreement is obtained at larger R values, as expected. 



4 Hybrid Mesons 



By hybrid meson, I mean a meson in which the gluonic degrees of freedom are 
excited non-trivially. I first discuss hybrid mesons with static heavy quarks 
where the description can be thought of as an excited colour string. I then 
summarise the situation concerning light quark hybrid mesons. 

Consider QQ states with static quarks in which the gluonic contribution 
may be excited. We classify the gluonic fields according to the symmetries of 
the system. This discussion is very similar to the description of electron wave 
functions in diatomic molecules. The symmetries are (i) rotation around the 
separation axis z with representations labelled by J z (ii) CP with represen- 
tations labelled by g(+) and u(-) and (iii) CJZ. Here C interchanges Q and 
Q, P is parity and 1Z is a rotation of 180° about the mid-point around the y 
axis. The CTZ operation is only relevant to classify states with J z = 0. The 
convention is to label states of J z = 0,1,2 by £, n, A respectively. The ground 
state (S + — g) will have J z = and CP = +. 

The exploration of the energy levels of other representations has a long 
history in lattice studies [TIJ. The first excited state is found to be the n n . 



This can be visualised as the symmetry of a string bowed out in the x direction 
minus the same deflection in the —x direction (plus another component of the 
two-dimensional representation with the transverse direction x replaced by y), 
corresponding to flux states from a lattice operator which is the difference of 
U-shaped paths from quark to antiquark of the form n — U. 

Recent lattice studies |16| have used an asymmetric space/time spacing 
which enables excited states to be determined comprehensively. These results 
confirm the finding that the U u excitation is the lowest lying and hence of most 
relevance to spectroscopy. 

From the potential corresponding to these excited gluonic states, one can 
determine the spectrum of hybrid quarkonia using the Schrodinger equation 
in the Born-Oppenheimer approximation. This approximation will be good if 
the heavy quarks move very little in the time it takes for the potential between 
them to become established. More quantitatively, we require that the potential 
energy of gluonic excitation is much larger than the typical energy of orbital or 
radial excitation. This is indeed the case |14|], especially for b quarks. Another 
nice feature of this approach is that the self energy of the static sources cancels 
in the energy difference between this hybrid state and the QQ states. Thus the 
lattice approach gives directly the excitation energy of each gluonic excitation. 

The U u symmetry state corresponds to excitations of the gluonic field in 
quarkonium called magnetic (with L PC = l + ~) and pseudo-electric (with 1~ + ) 



in contrast to the usual P-wave orbital excitation which has L PC = 1 . Thus 
we expect different quantum number assignments from those of the gluonic 
ground state. Indeed combining with the heavy quark spins, we get a degen- 
erate set of 8 states with J PC = 1~ , (T+ 1~+ 2 _+ and 1 ++ , + ~, l + ~, 2 +_ 
respectively. Note that of these, J PC = 1 h , + ~ and 2 + ~ are spin-exotic and 
hence will not mix with QQ states. They thus form a very attractive goal for 
experimental searches for hybrid mesons. 

The eightfold degeneracy of the static approach will be broken by various 
corrections. As an example, one of the eight degenerate hybrid states is a 
pseudoscalar with the heavy quarks in a spin triplet. This has the same overall 
quantum numbers as the S-wave QQ state (rjb) which, however, has the heavy 
quarks in a spin singlet. So any mixing between these states must be mediated 
by spin dependent interactions. These spin dependent interactions will be 
smaller for heavier quarks. It is of interest to establish the strength of these 
effects for b and c quarks. Another topic of interest is the splitting between the 
spin exotic hybrids which will come from the different energies of the magnetic 
and pseudo-electric gluonic excitations. 

One way to go beyond the static approach is to use the NRQCD approx- 
imation which then enables the spin dependent effects to be explored. One 
study [|l6j finds that the L PC = l + ~ and 1~ + excitations have no statistically 
significant splitting although the l +_ excitation does lie a little lighter. This 
would imply, after adding in heavy quark spin, that the J PC = 1 h hybrid was 
the lightest spin exotic. Also a relatively large spin splitting was found JT7| 
among the triplet states considering, however, only magnetic gluonic excita- 
tions. 

Confirmation of the ordering of the spin exotic states also comes from 
lattice studies with propagating quarks || [H| pl| which are able to measure 



masses for all 8 states. We discuss this evidence in more detail below. 

Within the quenched approximation, the lattice evidence for bb quarks 
points to a lightest hybrid spin exotic with J PC = 1~ + at an energy given by 
(m H -m 2S )r =1.8 (static potential pfl); 1.9 (static potential NRQCD 
2.0 (NRQCD [|HJ). These results can be summarised as (uih — rri2s) r o = 
1.9 ± 0.1. Using the experimental mass of the T(2S), this implies that the 
lightest spin exotic hybrid is at run = 10.73(7) GeV including a 10% scale er- 
ror. Above this energy there will be many more hybrid states, many of which 



will be spin exotic. A discussion of hybrid decay channels has been given p2 



The excited gluonic static potential has also been determined including sea 
quarks (Nf = 2 flavours) and no significant difference is seen [20]. Thus the 
quenched estimates given above are not superseded. 



I now focus on lattice results for hybrid mesons made from light quarks 
using fully relativistic propagating quarks. There will be no mixing with qq 
mesons for spin-exotic hybrid mesons and these are of special interest. The first 
study of this area was by the UKQCD Collaboration || who used operators 
motivated by the heavy quark studies referred to above to study all 8 J PC 
values coming from L PC = l + ~ and 1 h excitations. The resulting mass 
spectrum gives the J PC = 1~ + state as the lightest spin-exotic state. Taking 
account of the systematic scale errors in the lattice determination, a mass of 
2000(200) MeV is quoted for this hybrid meson with ss light quarks. Although 
not directly measured, the corresponding light quark hybrid meson would be 
expected to be around 120 MeV lighter. 

A second lattice group has also evaluated hybrid meson spectra with prop- 
agating quarks from quenched lattices. They obtain [|ll| masses of the 1~ + 
state with statistical and various systematic errors of 1970(90) (300) MeV, 
2170(80) (100) (100) MeV and 4390(80) (200) MeV for nn, ss and cc quarks 
respectively. For the + ~ spin-exotic state they have a noisier signal but ev- 
idence that it is heavier. They also explore mixing matrix elements between 
spin-exotic hybrid states and 4 quark operators. 

A first attempt has been made |19 to determine the hybrid meson spectrum 
using full QCD. The sea quarks used have several different masses and an 
extrapolation is made to the limit of physical sea quark masses, yielding a 
mass of 1.9(2) GeV for the lightest spin-exotic hybrid meson, which again is 
found to be the 1 h . In principle this calculation should take account of sea 
quark effects such as the mixing between such a hybrid meson and qqqq states 
such as T]7i, although it is possible that the sea quark masses used are not light 
enough to explore these features. 

The three independent lattice calculations of the light hybrid spectrum are 
in good agreement with each other. They imply that the natural energy range 
for spin-exotic hybrid mesons is around 1.9 GeV. The J PC = 1~ + state is found 
to be lightest. It is not easy to reconcile these lattice results with experimental 



indications |2T|] for resonances at 1.4 GeV and 1.6 GeV, especially the lower 
mass value. Mixing with qqqq states such as rjir is not included for realistic 
quark masses in the lattice calculations. This can be interpreted, dependent on 
one's viewpoint, as either that the lattice calculations are incomplete or as an 
indication that the experimental states may have an important meson-meson 
component in them. 



5 Conclusions 



Quenched lattice QCD is well understood and accurate predictions in the con- 
tinuum limit are increasingly becoming available. The lightest glueball is scalar 
with mass m(0 ++ ) = 1611(30)(160) MeV where the second error is an over- 
all scale error. The excited glueball spectrum is known too. The quenched 
approximation also gives information on quark- ant iquark scalar mesons and 
their mixing with glueballs. This determination of the mixing in the quenched 
approximation also sheds light on results for the spectrum directly in full QCD 
where the mixing will be enabled. There is also some lattice information on 
the hadronic decay amplitudes of glueballs. 

Evidence exists for a possible B S B S molecular state. 

For hybrid mesons, there will be no mixing with qq for spin-exotic states 
and these are the most useful predictions. The J PC = 1~ + state is expected 
at 10.73(7) GeV for b quarks, 2.0(2) GeV for s quarks and 1.9(2) GeV for 
u, d quarks. Mixing of spin-exotic hybrids with qqqq or equivalently with 
meson-meson is allowed and will modify the predictions from the quenched 
approximation. 
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